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Nonperturbative QCD Phenomenology and Light Quark Physics 

p. C. Tandy ^ 

'^Center for Nuclear Research, Department of Physics, Kent State University, Kent OH 44242, USA 

Recent progress in modeling QCD for hadron physics through truncated Dyson-Schwinger equations is reviewed. 
Special emphasis is put upon comparison of dressed quark propagators and the dressed quark-gluon vertex with 
lattice-QCD results. 
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1. Quark DSE solutions and Lattice Data 

The renormalised dressed-quark propagator, 
S{p), is the solution of the Dyson-Schwinger 
equation (DSE) 

S'\p) ^Z2ii-p + Zi m(/i) + S'(p, A), (1) 
wherein the dressed-quark self-energy is ^ 



g^D^,{k)^^^S{q)Tl{q,p).{2) 



E'(p,A) = ZiF 

Here Df^,j{k) is the renormalised dressed 
gluon propagator; k = p — q; T''^{q,p) is the 
renormalised dressed quark-gluon vertex; and 
ZiF(/i,A), Z2(/i,A) and Z4(/i,A) are, respec- 
tively, renormalisation constants for the quark- 
gluon vertex, quark wave function and current- 
quark mass. In Eq. @, -.^ d'^q/{2T:)'^ 
denotes an integral with a translationally in- 
variant ultraviolet regularization at mass-scale 
A. 

The general form of the dressed-quark propa- 
gator is S{p) = Z{p'^; fi^)/[i"f ■ p + M{p'^)] where 
Z{p^; and Af(p^) are the wave function renor- 
malization function and running mass-function 
respectively. The Z2 and Z4 are determined 
by solving the gap equation, Eq. J^l, subject 
to the renormalization condition ^"^(p)! 2- 2 — 
i"f ■ p + m[ii) at some large spacelike p.^. 

The gap equation, Eq. (P), is coupled to the 
DSEs satisfied by the gluon 2-point function and 



^We use a Euclidean metric where the scalar product 
of two four vectors is a ■ b = /^i—i '^i^i^ s-nd Hermitian 
Dirac-7 matrices that obey {7/j,7iy} = 2<5^^. 



the vertex 3-point function which in turn involve 
other n-point functions; a tractable model is de- 
fined by a truncation. At least one nonperturba- 
tive, chiral symmetry preserving truncation ex- 
ists |1I2| and the first term in that scheme is 
the renormalization-group-improved rainbow gap 
equation, wherein the self-energy, Eq. 0, as- 
sumes the form 



A" 



Q{k')Dl:'^(k)-^,S{q) -7., 



(3) 



where k — p — q and -D^''^''(fc) is the Landau gauge 
free-gluon propagator. 

In Eq. (PJ, t/(fc^) is an effective coupling, onto 
which has been mapped the combined effect of 
dressing both the gluon propagator and quark- 
gluon vertex. In the ultraviolet, G{k^) matches 
47ras~'°°P(fe^). It has been found that a one- 
parameter infrared form for S{k'^) can produce 
an excellent description for the ground state pseu- 
doscalar and vector mesons and their electroweak 
form factors and decays |3I4| . The advent of 
lattice-QCD data for the quark [S] and gluon |S] 
propagators provides an opportunity for an ex- 
amination of the QCD content of G{k'^)- 

A first effort has been made to unravel the 
infrared contribution to Gik^) from the dressed 
quark-gluon vertex [Jj. Since the lattice gluon 
propagator is infrared-suppressed, the required 
vertex will be infrared-enhanced to be empiri- 
cally successful. In that work, an effective ker- 
nel G{k^) is constructed from the quenched lat- 
tice Landau gauge gluon propagator 
and an effective infrared behavior of the vertex 
determined by requiring that the DSE reproduce 
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Figure 1. DSE solutions with vertex dressing fit to lattice data for the quark mass function for various 
current quark masses. 



the quenched lattice data ^ for the quark prop- 
agator. The kernel is given the form 



(4) 



where the dressed vertex is represented by 
7;.ri(fc2) with 

and T = — 1 > 1. The phenomcnological fac- 
tor v{k'^) is unity in the ultraviolet limit and the 
log factor combines with the leading log behavior 
of the gluon propagator to ensure that the entire 
kernel matches the 1-loop QCD behavior of the 
running coupling g{k'^) 47r^7m/ ln(fc2/A^cD). 
For the quenched, Landau-gauge study, Nj — 
0, du = 13/22 and 7„ = 12/(33 -2A^/) is the 
anomalous mass dimension. 

Studies of corrections to the rainbow trunca- 
tion show Fi to be the dominant amplitude of 
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the dressed vertex in respect to strength and ef- 
fect on observables. It is the only amplitude that 
is ultraviolet divergent at one-loop level. The 
current-quark mass, m{p), was fixed by match- 
ing the DSE and lattice results fS| for M{p'^) 
on > 1 GeV^ for each lattice set character- 
ized by ma. With a DSE renormalisation scale 
/i = 19 GeV and with m(/i) in GeV, we find 



.(6) 



In Fig. n we compare DSE solutions for M{p^) 
with lattice results. It is apparent that the lattice- 
gluon and lattice-quark propagators are easily 
correlated via a simple infrared vertex enhance- 
ment v{k'^) depicted in Fig. |21 Such vertex en- 
hancement has been anticipated for some time 
based on studies of the connection between DCSB 
and observables jHll()| . Note that v{Q) is finite; for 
the m{fj,) case considered in Sec. 3 in relation to 
lattice data for the vertex, v{k'^) = 13.0, where 
0.04 GeV^ is an effective infrared scale. 
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Figure 2. Vertex factor obtained in the chi- 
ral limit {solid curve) and with the three lowest 
current-quark masses from Eq. ©. v{0) is finite. 

With the dependence of the lattice data upon 
p and rn{iJ,) mapped into a DSE kernel, it is 
a vehicle for chiral extrapolaton. It is appar- 
ent from Fig. n that the linear extrapolation jJT] 
overestimates the chiral DSE mass function. We 
make a finer comparison in Fig. which dis- 
plays M(p^ — pjjj) versus m(/i) where pjj^ — 
0.38 GeV^, is the smallest value containing two 
lattice results for M{p'^) at the three lowest m{fi) 
values. The linear trajectory suggested by the 
lattice results provides a poor extrapolation to 

2 

m{ii) = 0, giving a result 40% too large . The 
linearly-extrapolated lattice data produces the es- 
timate 0: -(w)?GcV = (0.270 ± 0.027 GeV)3. 
However, the error is purely statistical. The 
systematic error, to which the linear extrapola- 
tion must contribute, was not estimated. The 
DSE-assisted chiral limit |Zj indicates that in 
quenched-QCD -(ot)?goV = (0.19 GeV)^, a fac- 
tor of two smaller than the empirical value (0.24± 
0.01 GeV)'^. Related to this, the quenched lattice 
data underestimates the chiral /° by 30% |7]. 

2. Beyond Ladder-rainbow 

Since the pseudoscalar states are dominated 
by chiral symmetry and its pattern of dynam- 
ical and explicit breakdown, most work on ex- 

^For an analysis in terms of the DSE chiral susceptibility 
see Ref. 1^1 . 
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Figure 3. Chiral extrapolation. Solid curve, the 
lattice-assisted DSE result; circles, lowest mass 
lattice data dashed-line, linear fit to the lat- 
tice data. 

tensions of the BSE kernel beyond the rainbow- 
ladder level have sought to preserve this symme- 
try. Particular extensions of the self-energy be- 
yond rainbow must be accompanied by related 
extensions of the ladder kernel so that the axial 
vector Ward-Takahashi is preserved. In the pres- 
ence of an empirical amount of DCSB, this iden- 
tity ensures the Goldstonc nature of the chiral 
limit pseudoscalars, and the physical mass pat- 
terns, independently of model details This 
may be implemented constructively if a model of 
the dressed quark-gluon vertex is defined in terms 
of Feynman diagrams. Several studies jll2ll4j 
have employed the Abelian-like series 

A^ 

ZiF—T^{k,p) = —-f,,- J Dpx{q)—-fp 

xS{k + q)—-f,S{p + q)—-fx + --- , (7) 

which has the ladder-rainbow truncation of 
Eq. (O as the first term. 

The corresponding Bethe-Salpeter kernel in 
this scheme is generated as the sum of terms pro- 
duced by cutting a quark line in the resulting 
diagrammatic representation of the quark self- 
energy. Due to the complexity of a BSE calcu- 
lation with a multi-loop kernel, most numerical 
implementations |1I2I14| have exploited the alge- 
braic structure that follows from use of the in- 




Figure 4. 1++ (ai//i) axialvector meson BSE 
eigenvalue from the ladder kernel [G = 0) and 
from the dressed vertex model (G > 0). 

frared dominant model D^^{q) oc GT^^{q)S*{q). 
For infrared-dominated physics this model effec- 
tively summarizes the behavior of more realis- 
tic models. A disadvantage is that it does not 
support finite quark-quark relative momentum in 
bound states. This is likely the reason why it does 
not bind the (P-wave) axial-vector states. 

Taking advantage of the fact that existing finite 
range ladder-rainbow models typically produce 
axial vector masses of about 0.9 GeV, Ref. ^B] 
has used the MN model only as an effective repre- 
sentation of the gluon exchange that implements 
vertex dressing via Eq. O . The residual effective 
gluon exchange is represented by an existing finite 
range model JSI- With a simple adjustment, the 
chiral symmetry-preserving BSE kernel is easily 
obtained. This does not implement the ultravio- 
let behavior of the QCD running coupling; it con- 
tributes typically 10% or less to meson masses 3 
and this level of precision is not our concern here. 

The eigenvalue behavior for the axial vector so- 
lutions in the 1"'""'" (ai//i) and 1"' [hi/ hi) chan- 
nels are displayed in Figs. 0] and |S1 In previous 
work, the ladder truncation, constrained by chiral 
data, is generally found to be 200-400 MeV too 
attractive for these P-wave states |17ll8ll6j . Our 
present results agree with this. The 1"'""'" channel 
shows a 30 MeV of attraction due to the effect of 
1-loop dressing added to the ladder kernel. How- 



Figure 5. 1"' {bi/hi) axialvector meson BSE 
eigenvalue from the ladder kernel (G = 0) and 
from the dressed vertex model (G > 0). 

ever, in the 1^ (bi/hi) channel, we find a re- 
pulsive effect of 290 MeV above the ladder kernel 
result, yielding a value close to m^^^. We are not 
able to compare these findings with previous work 
on vertex dressing since such P-wave states do 
not have solutions in the models considered pre- 
viously for that purpose |1I2I14| . Other studies 
of oi and bi based on the ladder-rainbow trunca- 
tion have used a separable approximation where 
the quark propagators are the phenomenological 
instruments jl9„20 21 , these studies find more ac- 
ceptable masses for both states in the vicinity of 
1.3 GeV. 

3. Nonperturbative Model of the Quark- 
gluon Vertex 

We denote the dressed-quark-gluon vertex for 
gluon momentum k and quark momentum p by 
igf^Trrip + k,p), where t'^ = X'^/2 and A'^ is an 
SU(3) color matrix. Through 0{g^), i.e., to 1- 
loop, the amplitude To- is given, in terms of Fig. 

by r,(p + k,p) = ZiF 7a + PA + + . . .. 

The color factors reveal two important 
considerations. The color factor of the 
(Abelian-like) term P^ would be given by 
t° t° = C'f = {N^ - l)/2Nc for the strong dressing 
of the photon-quark vertex, i.e., in the color sin- 
glet channel. The octet P^ is of opposite sign and 
is suppressed by a factor 1/{N^ ^1)- single gluon 




Figure 6. The dressing of the quark-gluon vertex 
at one loop: the AbeHan-Hke term F^, and the 
non-Abehan term F^^. 

exchange between a quark and antiquark has rel- 
atively weak repulsion in the color-octet channel, 
compared to strong attraction in the color-singlet 
channel. Net attraction for the gluon vertex (at 
least to this order) is provided by the non-Abelian 
term, which involves the three-gluon vertex: 
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Figure 7. The quark-gluon vertex amplitudes 
at zero gluon momentum and for quark current 
mass to(^ = 2 GeV) = 60 MeV. Quenched lattice 
data [H] is compared to the results of the DSE-Lat 
model 0. 

At fc = we note that the bare result 
^HAq^i)^ -d[D^\q)T^,iq)]/dq^ allows the 
total gluon content of the integrand of F^^ to be 
[dg'^ Do{q'^)/dqa]T^i,{q). The dressing provided 
by the combination F^ -I- F^"^ yields a vertex 
that satisfies the Slavnov- Taylor identity (STI) 
through 0{g^) |22|. 

Our nonperturbative model |22| for the dressed 
quark-gluon vertex is defined by extentions of the 
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two diagrams F^ -|- F^^ into dressed versions de- 
termined solely from an existing ladder-rainbow 
model DSE kernel that has 1-loop QCD renormal- 
ization group improvement. Two DSE models are 
employed. The first (DSE-Lat) 7 is the mapping 
of quenched lattice data for the gluon propagator 
into a continuum ladder-rainbow kernel with ef- 
fective gluon vertex as described in Sec. 1. In this 
sense, it represents quenched dynamics. The sec- 
ond (DSE-MT) ^ provides a good one-parameter 
fit to a wide variety of light quark meson physics; 
in this sense it represents unquenched dynamics. 
Both are implemented through the substitution 
Do{q^) G{q^)lq^ with bare quark propaga- 
tors replaced by rainbow DSE solutions. Note 
that this procedure makes a corresponding non- 
perturbative extension of T^J^^. Our justification 
is one of simplicity; no new parameters are intro- 
duced. 

The general nonperturbative vertex at 
fc = has a representation in terms of 
three invariant amplitudes; T„{p,p) — 
7m^i(p^) -4pM7-pA2(p2) -i2p^A3(p2); the 
lattice-QCD data ,8^ is provided in terms of 
the Ai(p^). A useful comparison is the Abelian 
Ward identity F^^(p,p) = -idS'\p)/dpa. 

In Fig. 13 we display the DSE-Lat model re- 
sults in a dimensionless form for comparison with 

3 

the (quenched) lattice data . The renormaliza- 
tion scale of the lattice data is fi — 2 GeV where 
Ai(/i) = 1, A{fj,) — 1. We compare to the lattice 
data set for which m(/i) = 60 MeV. The same 
renormalization scale and conditions have been 
implemented in the DSE calculations. Without 
parameter adjustment, the model reproduces the 
lattice data for Ai and A3 quite well. The Abelian 
Ansatz (Ward Identity) , while clearly inadequate 
for Ai below 1.5 GeV, reproduces A3. The DSE 
model A2 is evidently too weak, although the lat- 
tice data has large errors. The evident infrared 
enhancement at fc^^ 0.04 GeV^ is Xi{kf) ~ 2.3, 
a factor of 6 smaller than what the analysis of the 
quark propagator in Sec. 1 required. The non- 
Abelian term F^^ dominates to a greater extent 

^We note that in Ref. |S] both the lattice data, and the 
Abelian (Ward identity) Ansatz, for X3{p) are presented 
as positive. These two sign errors have been acknowl- 
edged 1^. 
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than what the ratio of color factors (—9) would 
suggest; it also distributes its infrared strength 
to favor Ai more so than does F^. The present 
model results could be summarized quite effec- 
tively by ignoring F^. Due to the definition of the 
two DSE models, their comparison |21] provides 
an estimate of the effects of the quenched approx- 
imation. The effects are comparable to present 
uncertainties in the lattice data. 

Recent work on a model of the quark-gluon ver- 
tex using a bare F^^^. vertex with DSE solutions 
for the gluon propagator and an Ansatz for dress- 
ing both internal quark-gluon vertices has pro- 
duced results similar to the present work, except 
that the m(/i) — 115 MeV case is considered |25| . 
Evidently the detailed infrared structure of F^^^. 
is not crucial to present considerations. 
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